Non-adiabatic pumping through interacting quantum dots 
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We study non-adiabatic two-parameter charge and spin pumping through a single-level quantum 
dot with Coulomb interaction. For the limit of weak tunnel coupling and in the regime of pumping 
frequencies up to the tunneling rates, Q < F/Zi, we perform an exact resummation of contributions 
of all orders in the pumping frequency. As striking non-adiabatic signatures, we find frequency- 
dependent phase shifts in the charge and spin currents, which opens the possibility to control charge 
and spin currents by tuning the pumping frequency. This includes the realization of an effective 
single-parameter pumping as well as pure spin without charge currents. 
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Introduction. Pumping is a transport mechanism which 
induces dc charge and spin currents in a nano-scale con- 
ductor in the absence of a bias voltage by means of a 
time-dependent control of some system parameters. Re- 
cently there have been several experimental works on 
pumping in nanostructures [H, 0, i, 0, i, 0, 0, Theo- 
retically, its interest lies in the possibility to investigate 
non-equilibrium phenomena induced by the explicit time- 
dependence of a nanoscale quantum system. A lot of in- 
terest has been devoted to the adiabatic regime, realized 
when the time-dependence of the parameters is slow in 
comparison to the characteristic time scales of the sys- 
tem, such as the dwell time. Many theoretical works 
have dealt with adiabatic pumping!: in systems with weak 
electron-electron interaction [§, 10, 11, 12, 1^, 3] as well 
as in systems where the Coulomb interaction cannot be 
treated in a mean- field approach [H, [H, [13, H El, 0, 
21, 22, 13, 13]. Pumping beyond the adiabatic limit, on 
the other hand, gives rise to larger pumping currents, 
facilitating the experimental investigation of this trans- 
port mechanism. Moreover, it adds another control pa- 
rameter, the pumping frequency Vt^ which can be used to 
steer the charge and spin currents. Pumping in the non- 
adiabatic regime is intrinsically a strong-non-equilibrium 
phenomenon and its theoretical description is challeng- 
ing. In the limit of weak Coulomb interaction, treated 
with a Hartree approach, a general theoretical framework 
can be based on the Floquet scattering matrix [25]. In 
systems with strong Coulomb interaction, no such gen- 
eral framework exists. As a paradigmatic system, we 
consider a single-level quantum dot. Pumping in this 
type of systems has been studied either in the adiabatic 

or in the opposite limit 
H, [29^. Typically, the 




regime 

of very large frequency [26', '2 
latter is studied in the context of photon-assisted tun- 
neling [lO, HH. On the contrary, in the present Letter 
we start from the low-frequency regime, including higher 
orders in the pumping frequency employing a diagram- 
matic real-time transport approach, which allows to in- 



clude Coulomb interaction and non-equilibrium effects. 
The transport quantities are then computed perturba- 
tively in the tunnel-coupling strengths. In the sequential 
tunneling regime all orders in the pumping frequency 
can be resummed. 

Model. We consider a single- level quantum dot, tunnel- 
coupled to two metallic leads. The Hamiltonian of the 
system is = i^dot + i^L + Hr + i^tunn- The dot is 
described by the Anderson impurity model 



(1) 



where ricr = d^^da- with the annihilation operator d^ for 
an electron with spin a =T,i. The eigenstates of i^dot 
are |x) with x ^ {0? corresponding to an empty, 

singly occupied with spin up or down, and doubly occu- 
pied dot respectively. The level position e{t) is periodi- 
cally modulated in time. Electrons in lead a = L^R are 
described by Ha = T^k,a ^c^kaclka^aka- The right lead 
can be ferromagnetic while the left one is non-magnetic. 
The spin polarization is pr = {ur^ - yR-)/{-UR+ + t^r-), 
where vrj^ i^R-) is the density of states at the Fermi 
energy for the majority (minority) spins in the right 
lead and < pr < 1. Dot and leads are coupled by 
^tunn = Ea'^al^jE/ca^La^^ + H.C.; the tunucl am- 
plitudes Va{t) vary in time. The tunnelling strength is 

^a{t) = Ea^c.a{t)/2 with Taa{t) = 2^^,, | (t) | ^ . The 

total tunnelling strength is T{t) = TL{t)-\-TR{t). No volt- 
age is applied: charge and spin currents arise only due 
to the periodic modulations of £{t) and Ta{t), denoted 
collectively by X{t) in the following, with frequency ft. 
Non-adiabatic pumping. The dot is described by its re- 
duced density matrix pdot(^)- In the present case, the 
dynamics of the diagonal and off-diagonal elements of 
Pdot(^) are decoupled, i.e., we can restrict ourselves to 
study the occupation probabilities, Px{t) = (xIPdotlx)? 
whose time evolution is governed by a generalized mas- 
ter equation 



'it) = f 

J — C 



dt'V^{t,t')-'P{t'), 



(2) 
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where P = {Pq, , Pi, Pd)^ • The kernel W{t,t') func- 
tionahy depends on X{t). We perform a series expan- 
sion in powers of the pumping frequency. Generahz- 
ing the adiabatic expansion of Ref. [2l|, we keep all or- 
ders in Q. For this, we write P(t) = ^ 



k>0 ^ t 



and 



= E/c>o - where the superscript (k) 

indicates terms of order ft^. In W^^^ {t — t'), the time de- 
pendence of ah parameters is expanded to order k around 
the final time t and only terms of order k in the time 
derivatives are retained. Expanding P(t') in Taylor se- 
ries and performing a Laplace transform of the r.h.s. of 
Eq. ([2]), we obtain 



n=0 p,(7,/c=0 



m 



(3) 



where S^W^^^ 



= lim^^o+ d^Wr {z)/dz^ and 
W^^^(z) = n/^^e-^(^-^')W^^^\t - tO^it' is the Laplace 
transformed kernel. How to count the time derivatives of 
P in this expansion depends on the considered frequency 
regime. In this paper, we consider the regime hTt < F, 
for which the system quickly relaxes to an oscillatory 
steady state with the frequency of i.e., each time 

derivative introduces one power in Vt. 

In addition to the expansion in frequency, we perform a 



systematic expansion of W^^^ = Xlj>i and PJ 

X]j>-fe -^tj' powers of the tunnel-coupling strength F. 
The order in F is indicated by the subscript j. Upon 
substitution into Eq. ([3|), the orders of and F on both 
sides are matched, giving rise to a hierarchy of coupled 
equations for V).J . This matching requires the expansion 

for P^^^ to start from F~^. In the rest of the paper, we 
concentrate on the limit of weak tunnel coupling, i.e., we 
compute the kernel to first order in F. In this case, the 
hierarchy of equations reduces to 



o = w^^pl^npl^. 



w 



(0) 



^,-(/c+l) • 



(4) 



for /c > 0. Remarkably, only the instantaneous kernel 
w|^^\ corresponding to freezing the time evolution of 
X{t) at time t, is needed. The rules for evaluating this 
kernel are given in Ref. [24]. We can solve for P[^2/. re- 
cursively starting from the instantaneous term p[^q . The 
charge and and spin currents in the the left lead can be 



expanded in powers of Q as well, /^(t) 



= (5,5'). The k-th contribution is given by I^^\t) = 



^^,L(o) p(fc)^^ ^-^^ =e, cs = h/2, and 

= (1, . . . , 1). The symbols wf'/"^^^ stand for the cur- 
rent rates, which take into account the number of elec- 



trons transferred to the left lead [2J. In steady state, the 
average pumped charge and spin currents per period are 



chosen for calculating the average pumped currents per 
period, since the current continuity equation is fulfilled 
for each order of the expansions. 

The procedure outlined above is very general, pro- 
vided the weak-coupling limit. To proceed in our spe- 
cific model, we introduce the dot charge (in units of 
e) NQ{t) = (0,1,1,2) • P{t) and spin (in units of h/2) 
Ns{t) = (0,1,-1,0) • P(t) and their deviations from 
the instantaneous values AN^ = N^{t) — N^^\t) with 

Ng\t) = 2/[s(i)]/{l+/[s(t)]-/[s(t)+C/]}andiVf (t) = 
0. A resummation of Eqs. dH) yields 

AN^+T-\t)AN^+PnT^\thR{t)AN^ = -N^°\t), (5) 

where ja{t) = Ta{t)/T{t) and ^ = 5 if ^ = Q and vice 
versa. The two time scales r^(t), defined by 



hTsHt) 



r{t){l + f[e{t)]-f[e{t) + U]} (6) 
r{t){l-f[s{t)]+f[s{t) + U]}, (7) 



mare the instantaneous charge and spin relaxation times 
for = 0, with f{E) being the Fermi function. Finally, 
the charge and spin currents can be recast in the form 



i^{t) = -c^T^HthLmm) ■ 



(8) 



Equations ([5]) and (|8|) constitute the main result of 
this Letter. They allow to evaluate the non-adiabatic 
pumped charge and spin currents for frequencies hft <r. 
Solutions for both = and Pr are discussed in 
the following. As a specific pumping model, we choose 
^L{t) = fL + AFLsin ((^t) and6(t) = £^Aesm{nt ^ 
with (/) the pumping phase. We take F^^ as time inde- 
pendent. We focus chiefiy on the case of weak pumping 
ATl, Ae <C F which allows a analytical treatment. How- 
ever, we want to stress that Eqs. (jSHH]) are not restricted 
to this regime. Numerical results for strong pumping will 
be discussed in the last part of this work. 

In the weak-pumping regime, the charge and spin cur- 
rents have the form 



= Jf"^sin((/) + A(/)^), 



(9) 



where both the amplitude X^^^ and the phase shift Acf)^ 
are odd functions of (7, i.e., when expanding the currents 
in powers of the pumping frequency, all the odd powers 
of the frequency Q are proportional to sin(0) while the 
even powers are proportional to cos(^). The analytical 
expressions for the ^l-dependent amplitudes are, for ar- 
bitrary value of the spin polarization pR, rather lengthy 
and we do not report them here. 

The zeroth order in Q describes the instantaneous equi- 
librium, for which both the average charge and spin 
current vanish. Adiabatic pumping corresponds to ex- 
panding the currents to first order in which leads to 
{dX?^^^/dft)\Q=oftsin(l). The non-adiabatic con- 



rad 



= /o^ dtmt). It is immaterial which barrier is tributions to the pumped charge and spin do not only 



3 



introduce an 11-dependence of the amplitude, but change 
the pumping behavior qualitatively. First, phase shifts 
are generated, i.e., pumping is also possible for 
^ = or TT, which corresponds to single-parameter pump- 
ing. Second, the phase shifts for charge and spin pump- 
ing may differ from each other. As we will see below, 
this allows pure spin pumping [32] without charge pump- 
ing. Phase shifts are a general feature of non-adiabatic, 
driven quantum dynamics and are found in disparate con- 
texts, from circuit-QED [33] to driven optical lattices 
and stochastic quantum resonance [35|]. Molecular sys- 
tems [sg} and nano-electomechanical systems [s^ are es- 
pecially interesting in connection to our technique. 
Non-magnetic case. In the limit of non-magnetic leads, 
PR = 0^ the expressions for the charge and spin currents 
simplify substantially. Then, the spin currents vanishes, 
^5 = 0, and for the charge current we find 



A0Q 



ArLA£G""!^fQ[l 
— arctan (f^fg) . 



(10) 
(11) 



Here, we used the expression for the linear DC con- 
ductance r^li" - fj ^LTB f(s)\l-f(s+Um-f(s)+f(s+U)] 

auctance u _ ^ p - , 
which, for ksT <C U has two maxima at £ = ksT \ii2 
and e = —U — ksT \n2. 




FIG. 1: (Color online) (a) Charge- current amplitude Iq"^"^ 
(thin solid), its adiabatic approximation Xg^^' (thin 
dashed), and phase shift A(/)g (thick) as a function of hQ/T. 
(b) Pumped current Xq as a function of hQ/T, for = 0.27r 
(solid), = 0.47V (dashed), = O.Gtt (dash-dotted), (c) 
Pumped charge current 1q as a function of 0, for hQ/F = 0.1 
(solid), 0.5 (dashed) and 0.8 (dash-dotted). All currents 
in units of e[3/W l^s/H. Other parameters: e — /csTln2, 
U 9f , kBT 3f , PR 0, and 7^ 1/2. 



In the following, we concentrate on the peak at £ = 
kBThi2. In Fig.flfa) the charge-current amplitude Tq^^ 
and the phase shift A0q are shown (solid lines) as a func- 
tion of hQ./V. The amplitude is reduced as compared to 
the adiabatic approximation (dashed line). In Fig. [TJb) 
the pumped current Xq is plotted as a function of the 
frequency Vt for different values of the pumping phase 
(j). Deviations from the adiabatic approximation, simply 
given by the tangents at 1] = 0, become most pronounced 
as the pumping phase (j) approaches zero: in this limit the 
phase shift A(/)q, a striking signature of non-adiabatic 
pumping, is most important. Figure [Tfc) shows Xq as a 
function of the phase (j) for different values of hTt/V. The 
growing relevance of the phase shift A0q as the pumping 
frequency is increased is apparent. As commented above, 
its most important consequence is the possibility to ob- 
tain an effective one-parameter pumping for ^ = 0, tt. 

The frequency dependence of the phase shift A^q is 
that of a low-pass filter with cutoff frequency fg^. It sig- 
nals the reduced ability of fast pumps to transfer charge 
due to the non-zero response time of the dot to charge 
excitations, tq. 



0.007 




-0.006 



hQ./T 



FIG. 2: (Color online) (a) Spin-current amplitude X™'''' 
(thick) and its adiabatic approximation x™^^'^^ (thin) as a 
function of hQ./V for different values of pR. (b) Phase shift 
A(/)5 as a function of hQ./V for different values of pr.{c) 
Charge and spin currents Xq^s as a function of HQ/T with 
PR = 0.4 and = 7r/9 for weak pumping (thick) and strong 
pumping (thin, with AFl = Ae 0.4f). Spin currents 
in units of ;^Ari:A£/2. Other parameters: s = /csTlog2, 
U = 9f , ksT = 3f , and ^r = 1/2. 



Ferromagnetic case. In the weak-pumping regime for 
Pi? 7^ 0, we resort to a graphical presentation of the av- 
erage pumped currents. In Fig. [2l^a) and (b) the spin- 
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current amplitude 2"^^^ and phase shift A^^ is shown 
as a function of the frequency Q. Both exhibit a strong 
dependence on the polarization pji of the ferromagnet. 
By contrast, the charge amplitude and phase shift (both 
not shown here) are almost insensitive to pr 0. In 
particular, A(j)Q is essentially given by Eq.fTTland charge 
transfer is still dominated by the cutoff frequency fg^. 
On the other hand, the intricate, coupled dynamics of 
spin and charge results in the spin phase being given 
by A(j)s = arctan(l]/21]o) — Xlg=± ^^ctan(l]/21]g). For 
= 1/2, the cut-off frequencies are = p\fg^ and 
^± = Tq' + ± ^{TQ'-Ts'y+pUfQfs)-'. The 
competition between and the "hybridized" Q± gives 
rise to the non-monotonic behavior as well as the sign 
change of A^^. The most striking consequence of non- 
adiabaticity in the ferromagnetic case is the possibility 
to obtain a pure spin current, Xs with = 0, by 
tuning Q such that A0q(1]) = 0. This behavior is il- 
lustrated in Fig. [2jc), thick lines. The adiabatic approx- 
imation, missing the phase shifts, is unable to capture 
this effect. Pure spin pumping, as all the effects induced 
by the ferromagnetic lead, are enhanced for ^ 1. 
We conclude discussing the experimental observability of 
the effects discussed above. Currents in the weak pump- 
ing regime are generally small. However, our theory is 
valid also for strong pumping where currents are larger 
and we have checked that all the results are not quali- 
tatively modified. As an example. Fig. [2](c) (thin lines) 
shows the charge and spin currents for ATl = As = 
OAT. For a temperature T of the order of 1 K {ksT ^ 
90/i eV) and f = ksT/S , one has e(3ATLAe/h ^ 500 
pA as the current unit. Clearly, pure spin pumping is 
still present. The charge current, ranging between -3 pA 
and 1 pA, is easily measurable. However, care needs to 
be exerted to minimize rectification effects due to stray 
capacitances and to distinguish them from pumping. 
Conclusions. We have studied non-adiabatic two- 
parameter charge and spin pumping through an inter- 
acting quantum dot, tunnel coupled to metallic leads. In 
the sequential-tunneling regime, for frequencies smaller 
than the tunnel rates, all orders in the pumping frequency 
can be resummed exactly. We find frequency-dependent 
phase shifts for the pumped charge and spin current, ne- 
glected in the adiabatic approximation. They allow an 
effective single-parameter pumping, i.e., finite pumping 
currents when the two pumping parameters oscillate in 
phase or anti-phase. When one lead is ferromagnetic, the 
pumping frequency turns out to be an important param- 
eter to control the spin current. In particular, due to 
the non-adiabaticity of pumping, pure spin currents can 
be generated. The effects discussed above can be exper- 
imentally verified in the strong-pumping regime, where 
the pumped currents are sizable. 
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